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Abstract. 

We study the phenomenology of supersymmetric models that explain neutrino masses 
through the spontaneous breaking of R-parity, finding strong correlations between the decays of 
the lightest neutralino and the neutrino mixing angles. In addition, the existence of a Goldstone 
boson, usually called Majoron (J), completely modifies the phenomenology with respect to 
the standard picture, inducing large invisible branching ratios and charged lepton decays, like 
jj, — > eJ, interesting signals that can be used to constrain the model. 



1. Introduction 

Nowadays it is well stablished that neutrinos have non-zero masses and mixing angles. In fact, 
oscillation experiments [H EJ [Sj HI [5] have become more precise in the last years, allowing us to 
measure the involved parameters with good accuracy. Moreover, global fits [6| have shown the 
existence of two large mixing angles, 9i2 and 023, one small (perhaps zero) mixing angle, ^13, 
and two different mass scales, given by the squared mass differences Am^2 "^"^23 • 

Therefore, since the minimal models do not include neutrino masses, we must enlarge them in 
order to explain this data. Many proposals can be found in the literature with that purpose. The 
most popular ones are those based on the different variations of the so-called seesaw mechanism 
El [Hi [TOl [11] , in which the existence of a high energy scale can naturally explain the smallness 
of neutrino masses. However, it is impossible to test these models directly, and only with 
some additional assumptions one might be able to get indirect insights thanks to forthcoming 
experiments [HI |13l [HI [IS [l6] . 

Consequently, it is worth to consider alternative models, in which neutrino masses are 
generated at the electroweak scale. In this case one can also find many different ideas. One 
possibility is to get neutrino masses through radiative corrections, like in the Zee model [T7] or 
in the Babu-Zee model |18j . two good examples of models with a rich phenomenology at current 
colliders like the LHC. On the other hand, in the context of supersymmetry, the violation of 
R-parity {]^p ) [19] is an interesting alternative, since it can provide connections between the 
decay of the lightest supersymmetric particle, the LSP, and neutrino physics [20l [2111221 [23 l [M] . 

There are many variations of J^p models. The simplest versions break R-parity including 
explicit terms in the superpotential that break either lepton or baryon number. Following that 
idea, of particular interest for neutrino physics are bilinear R-parity violating models (BRpV) 
|25j . in which lepton number violating terms eHL are introduced. The mixing between the 



neutrinos and the rest of neutral fermions induce small neutrino masses if the dimensionful e 
parameters are small compared to the SUSY scale, giving a TeV scale version of the seesaw 
mechanism. 

However, there is no explanation for the smallness of the e parameters, since they are SUSY 
conserving terms, and therefore their natural scale is the Planck scale. Spontaneous violation 
of R-parity {s-]^„ ) [261 [27] can give a natural solution, generating the bilinear terms at the 
electroweak scalq^ while keeping the nice features of the simplest models. Moreover, in addition 
to the usual signatures related to LSP decays, the spontaneous breaking of lepton number 
implies the existence of a Goldstone boson, the Majoron (J), which can strongly modify the 
phenomenology. For example, as it was already shown in |29[ I3U|. the invisible decay channel 
Jv can be dominant, leading to a possible confusion with the usual conserved R-parity 
models. 

In this paper we will focus on the phenomenology associated with the decays of the lightest 
neutral fermion, and how we can extract information connected to neutrino physics, for different 
scenarios: (i) a bino-like LSP, very common in the literature for mSUGRA models, and (ii) a 
singlino-like LSP. In particular we will confirm the result in [291 130| . where the difficulties to 
distinguish the model from R-parity conserved in special regions of parameter space are shown. 
In connection with that, exotic muon decays like /x — > e J and ^ — > e J7 are enhanced in the same 
regions of parameter space where the invisible decays of the LSP are dominant, thus providing 
information complementary to accelerator searches. 

The paper is organized as follows: in Section [2] we present the model with special emphasis 
on neutrino physics. In Section [3] we discuss production and decays of the lightest neutralino at 
the LHC. In particular, we will stress the importance of the different correlations between ratios 
of neutralino decay branching ratios and neutrino mixing angles. In Section H] we study exotic 
muon deays involving Majorons and how to use them to constrain the model. Finally, we draw 
our conclusions in Section [5l 

2. The model 

2.1. Model basics 

The model we consider |27) contains three additional singlet superfields, namely, D^, S and <I>, 
with lepton number assignments of L = —1,1,0 respectively. Then, taking into account two 
basic guidelines: (i) R-parity is conserved at the level of the superpotential, and (ii) We do not 
include terms with dimensions of mass, we can write down the following superpotential: 

+ li',Li^'Hu-hoHdHu^ + M>d'S+^^^ (1) 

It must be noticed that the introduction of the $ superfield addresses the ^u-problem of 
supersymmetry a la NMSSM [31jj. Similarly, the right-handed neutrino superfield v'^ will play a 
similar role in the spontaneous generation of bilinear terms, through its Yukawa coupling with 
the left-handed neutrino superfield. For both the inclusion of the second guideline is needed. 
Finally, the soft supersymmetry breaking terms of this model can be found in (32) . 

When electroweak symmetry breaking occurs, several scalar fields acquire vevs [33]. Besides 
the usual Higgs bosons 

^ Recently, other models in which the bilinear terms are generated spontaneously have appeared in the literature. 
As a simple but interesting example, the so-called ^i^SSM [28] has many things in common with s-I^p , promising 
a rich phenomenology at colliders. 



{H'd) = ^ iH'u) = ^ (2) 

we also have 

<*>-5 <->-^ <^>-^ <^->-^ 

As a consequence of this vacuum structure, an effective // term appears, fi = ^o^j which is 
naturally at the electroweak scale. Moreover, R-parity is spontaneously broken by the vevs vr, 
vs and vl-, and bilinear terms = generated as well. In this model their smallness 

is naturally explained by the smallness of the Yukawa couplings of the neutrinos, /ij,, which are 
also connected to the vevs of the left-handed sneutrinos, vl^, through the tadpole equations [2Tj. 
Therefore, we will consider in the following that both and vl- are small parameters compared 
to the usual SUSY scale. 

This spontaneous breaking of lepton number leads to the appearance of a Goldstone boson, 
usually called Majoron, J. The existence of a massless Majoron is usually seen as a bad feature 
of models, since the first s-^p models [26] were ruled out by LEP and astrophysical data [Ml |35] 
because of the doublet nature of this Goldstone boson. However, in this model R-parity is 
broken by the vevs of singlets, and therefore the Majoron has singlet nature, avoiding the strong 
constraints. In fact, it is possible to find a simple approximation for the Majoron 

J.(0.0,f,0,^f,-^5) (4) 



where V = \Jv'^ + v'g and terms of order ^ have been neglected. As we will see, the presence 

of this massless particle can strongly modify the phenomenology at colliders from the standard 
picture, introducing invisible decay channels which can have large branching ratios. Moreover, 
it can also have an impact on charged lepton decays, opening the possibility to processes like 
li — > IjJ and li IjJ^. 

2.2. Neutrino masses 

Once lepton number is broken neutrinos get masses through neutralino-neutrino mixing. In the 
basis (ip^)'^ = {B^, W^, Hd, Hu, fe, t'/x, i^t, i^'^, S, $) the 10 x 10 neutral fermion mass matrix can 
be written as [32l [36] 

/ Mh m3x7 

Mn = I (5) 

V ^Lr 

The sub-block Mh is the mass matrix of the seven heavy states while 1113x7 is the matrix 
that mixes them with the light neutrinos. Their exact expressions can be found in [30]. 

Since neutrino masses are much smaller than all other fermion mass terms, one can find the 
effective neutrino mass matrix in seesaw approximation 

and, after some straightforward algebra, can be cast into a very simple form 

- i'^tZ)ij = O'^i^i + K^i^j + ^3^i) + ceiej (7) 



mr.V. = -m3x7 • Mh ^ • m|'x7 (6) 



where Aj = eiVd + ^^vl^ are the so-called alignment parameters, that also appear in BRpV, 
and a, b and c are different combinations of the parameters of the model. 



It is numerically found that the 6-term is smaller than the other two, leaving us with two 
pieces to fit neutrino data. Then, two important conclusions can be learnt from equation ([7]): 

• The matrix has two non-zero eigenvalues eigenvalues at tree-level, as opposed to BRpV, 
where radiative corrections are needed to generate a second mass scale. 

• There are two possibilities to fit neutrino data: 

— Case (cl): A generates the atmospheric scale, e the solar scale 

— Case (c2): e generates the atmospheric scale, A the solar scale 

Since there are no a priori reasons for choosing one of these two options, we will consider 
both in the following. In any case, the l^p parameters and Aj have to be small in order to 
predict a correct absolut scale of neutrino mass. For typical SUSY masses order 0(100 GeV), 
\A\/iJ.'^ ~ 10~^-10~^. If some of the singlet fields are light, i.e. have masses in the range of 
0(0.1 - few) TeV, also \ei/iJ.\ can be as smaU as |el/^ ~ lO-^-lQ-^ 

3. Neutralino phenomenology 

In this section we discuss the phenomenology of a neutralino LSP in s-I^p at future colliders, 
focusing on the most important qualitative features. All numerical results shown below have 
been obtained using the program package SPheno [37], extended to include the new singlet 
superfields of this model. For details about the numerical procedure see |30j . 

3.1. Bino vs Singlino 

In principle, it is impossible to give a definite prediction for the nature of the LSP. This is due 
to the many unknown parameters in the spectra of the seven heavy neutral fermions. But there 
are general properties, present in most parts of parameter space: 

• There are four states very close to the MSSM neutralinos. In mSUGRA models the lightest 
of these four states is the Bino, with a mass given by Mi. 

• The states u'^ and S form a quasi-Dirac pair, the so-called Singlino, Si^2 — '^(^^ "F 'S')- 

• The remaining state is the phino, also present in the NMSSM. 

We are going to focus on a neutralino LSP, studying two different cases concerning its nature, 
(i) A Bino-like Xi^ ^ ^ typical mSUGRA scenario, and (ii) a Singlino- like Xi^ ^ novel but 
natural possibility. The parameter that determines in which case we are is the u'^ — S element in 
the neutral fermion mass matrix. If it is smaller than the Bino mass the Singlino will be lighter 
than the rest of neutralinos. Taking the expressions for the mass matrices (see reference [30] for 
the formulas) one finds that this condition reads M^^cg = -^hv^ < Mi, what can be naturally 
achieved. 

3.2. Production 

Since neutrino physics requires that the I^p parameters are small, the production cross sections 
are very similar to the corresponding MSSM values [38]. Therefore, we expect standard decay 
chains, with gluinos and squarks directly produced at the LHC and a neutralino LSP appearing 
as part of the final state (actually, it will further decay thanks to the breaking of R-parity). For 
example, a decay chain like 

q + B ^ q + Si^2 + J ^ ■■■ (8) 

will produce singlinos as an intermediate state, before starting with the ]^p decays. 

In conclusion, the production of the lightest neutralino is guaranteed, both for a Bino-like 
LSP and for a Singlino-like LSP. Moreover, we see that singlinos can be produced at accelerators, 
although their direct production cross section is negligible, allowing us to study their decays. 



3.3. Decays 

With broken R-parity the hghtest super symmetric particle decays. In figure [T] we show the 
most important decay channels for a neutralino with m^o > my^± . Note the transition between 

Singhno LSP (left) and Bino LSP (right) for ^hv^ ~ Mi (Mi ~ 98 GeV in SPSla' standard 
point) and how the different branching ratios are modified depending of the nature of the lightest 
neutralino. 

For low values of the LSP mass the channel J / Sj + v is usually the most important. The state 
Sj is a rather light singlet scalar, called the scalar partner of the Majoron in [36], that decays to 
a pair of Majorons with a branching ratio close to 100%. Therefore, this channel represents 
invisible final states, since the Majoron escapes detection, and it has a sizeable branching 
ratio even for a relatively high vr. We will come back to this important point later. Next 
in importance are the final states involving and charged leptons. The relative size of the 
branching ratios for the final states We, Wfi and Wt depends on both, (a) the nature of the 
LSP and (b) the fit to the neutrino data. Finally, the model predicts 

E,i3r(xg^Z° + r/,) _ g 
2J2iBr{x^^^W+ + 1^) 4cosHw 

with g being a phase space correction factor, with — > 1 in the limit m^o 
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Figure 1. Branching ratios for the most important decay modes of the lightest neutralino 
versus '^hv^ for a specific, but typical example point, in which the MSSM parameters have 
been adjusted such that the sparticle spectrum of the standard point SPSla' is approximately 
reproduced. The singlet parameters have been chosen randomly, vr = vs = 1 TeV, e and A 
have been fitted to neutrino data, such that A generates the atmospheric scale and e the solar 
scale. 

For the case of m^o < m^^^^i , one finds that invisible decays also play a very important role, 
with large branching ratios, which can be even dominant in some parts of parameter space. The 
rest of decay channels are three body decays mediated by virtual gauge bosons. See [HU] for the 
details. 

The induced invisible decays are the most important consequence of the existence of the 
Majoron. In the previous examples there are regions of parameter space where they are 
dominant, but still allowing for a sizeable branching ratio to visible final states. 
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Figure 2. Sum over all at least partially visible decay modes of the lightest neutralino versus 
vji in GeV, for a set of values f $ = 10 — 40 TeV for the mSUGRA parameter point mo = 280 
GeV, mi/2 = 250 GeV, tan/3 = 10, = -500 GeV and sgn{^i) = +. To the left Xi ^ B; to 
the right Xi — <S- 

Figure [2] shows the sum over all at least partially visible decay modes of the lightest neutralino 
versus vr in GeV, for a set of v$ values = W — 40 TeV for the mSUGRA parameter point 
(mo = 280 GeV, m^/a = 250 GeV, tan^ = 10, Aq = -500 GeV and sgn(//) = +). This point was 
constructed to produce formally a Q^oh^ ~ 1 in case of conserved R-parity, much larger than the 
observed relic DM density [30]. The left plot shows the case Xi — B, the right plot Xi ~ <S- The 
plot demonstrates that the branching ratio into B ^ Jv does depend strongly on the value of vr 
and to a minor extent on Vi^. Lowering vr one can get branching ratios for the invisible decay of 
the B very close to 100%, thus a very MSSM-like phenomenology, since the Majorons mimic the 
usual missing Et signal of the MSSM. This feature is independent of the mSugra parameters, 
see the correspoding figure in ^29j. In this case large statistics becomes necessary to find the 
rare visible neutralino decays, which prove that R-parity is broken. The inconsistency between 
the calculated Q-^oh? and the measured VlcDMh'^ might give a first indication for a non-standard 
SUSY model. Also, as we will comment in section HI exotic muon decays are enhanced for low 
values of vr, and therefore they can provide additional tools to test the model in that region. 

Figure [2] to the right shows that the case Xi — has a very different dependence on vr. We 
have checked that this feature is independent of the mSugra point. For other choices of mSugra 
parameters larger branching ratios for Br (5 — > J + z/) can be obtained, but contrary to the Bino 
LSP case, the sum over the invisible decay branching ratios never approaches 100%. 

Another important property of the decays of the lightest neutralino is their strong correlation 
with neutrino physics. Since LSP decay and neutrino masses are consequences of the same 
physical input, the breaking of R-parity, some connection between them is expected. In 
particular, one can use the decays to Wli (with the W boson real or virtual) to measure neutrino 
mixing angles. 

This is possible due to the good approximations that can be made in the — — 
couplings, which are found from the general expressions for the x^ — — x^ vertices (see |30j 
for explicit formulas). Using these approximations one can show that 

Br{B^Wli) (X (10) 
Br{S^Wk) (X ei^ (11) 




where the proportionality constant is independent of the generation index 'i' for both cases. 



This feature is shown in figure [3] 
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Figure 3. Ratios of branching ratios Xi ~^ Wli and their correlations with I^p parameters. 
To the left: Ratio _^^|^o~^|^|^^^ versus (Ae/A^)^ for a Bino LSP. To the right: Ratio ^^|^oZ^!^|^| 
versus (e^/e,-)^ for a Singlino LSP. All points with niLSP > mw. 



Since the structure of is given by A and e, this implies correlations between some 

combinations of branching ratios and neutrino mixing angles. 

The correlations depend on the nature of the LSP, Bino or Singlino, and the fit to neutrino 



data, (cl) or (c2). In figure[l]we give a representative example, showing the ratio 



as 



a function of tan^ OAtm for a Bino LSP. Using the current experimental range for the atmospheric 
angle one can predict an allowed range for the ratio, which has to be within [0.4, 2.1] to be fully 
consistent with the model. 
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Figure 4. Ratio Br{xi^-rW) ^^^^^^ ^^^^ QAtm for a Bino LSP. Vertical lines are the So" c.l. 
allowed experimental ranges, horizontal lines the resulting predictions for the fit (cl). 

For other cases see [30]. It must be stressed that for the case m^o < 'm-^± the correlations 
are still present, with virtual W bosons instead of real ones, and similar results can be found 
with the channels x? ~^ hll' ■ 

In conclusion, these correlations make the model testable at the inminent LHC. However, since 
we don't know whether case (cl) or case (c2) is realized, the decay of the lightest neutralino is 



not sufficient to determine the nature of the LSP. We need to reconstruct the complete decay 
chains and use kinematical variables to obtain some information about the intermediate states. 

4. Exotic muon decays 

As it has been shown in section [3l low values of vr enhance invisible Bino decays, and thus makes 
it difficult to distinguish the model from conserved R-parity. If Nature has chosen that region 
of parameter space we will need really high statistics to prove R-parity breaking at colliders. 
For that reason, it is interesting to have other observables to test the possible existence of the 
Majoron. 

Exotic muon decays involving Majorons [H], like — > e J and /u eJ^, might provide such 
experimental input, since they are also enhanced for low values of vji. This relation can be 
directly seen from the coupling li — J — Ij, which has the form 

O'H^j = {OZjY (12) 

and the C coefficients are different combinations of MSSM parameters. Note that the coupling 
is divided hy V ^ vr, and thus it is larger for low values of vr. 

By measuring deviations from the standard decays of the muon one can discover these decay 
channels or, at least, put bounds on the parameters of the model. In particular, if these processes 
are not observed at low energy experiments, their results will be translated into lower bounds 
for vr. If those bounds are in conflict with the non-observation of visible decays of the lightest 
neutralino the model will be in trouble. 

Actually, there are current bounds on the branching ratios of the decays fi ^ eJ and fj, — > e J7. 
On the one hand we have results by experiments studying muon decay and looking for fi ^ e'j, 
which indirectly can provide information about these signals. As it was already studied in [12], a 
visible peak over the background signal ^ — > eui? would mean the discovery of ^ — > eJ. However, 
the current bound Br{fi eJ) < 2, 6 • 10^^ , taken from that reference, does not directly apply 
to our model, since the authors consider an isotropic familon emission, which is not the case 
for a Majoron with parity violating couplings. On the other hand, astrophysics puts a worse 
but more reliable bound thanks to the possible influence of the Majoron in the cooling of red 
giants. In [HI] it is estimated that the coupling e - e - J has to obey Qeej < 10^^^ in order 
to leave important stellar properties unchanged. This can be translated into an approximate 
bound Br{fi eJ) < (few) ■ 10~^. Finally, for the case of /u — > eJ^, the reference [43] gives the 
bound Br{^ — > eJ^) < 1.3 • 10^^, although is is restricted to a small part of the phase space. 

In figure [5] the branching ratio of — > e J is shown as a function of vr for different values of 
v^. As previously explained, for low values of vr the branching ratio can be extremely large, 
with values clearly above the current bounds. In addition, there is an indirect dependence on 
f$ through neutrino masses. The parameters o, b and c in equation d?]) are dependent on f$ 
and they are strongly decreased for high values of this parameter. Then, if f$ is large one needs 
also high values for the parameters and Aj in order to correctly fit neutrino masses, thus 
indirectly increasing the coupling ()12p . 

The decay /i — > e J7 can be also of interest for our purposes. Many of the current experiments 
measuring muon decay are designed to look for photons in the final state, since /i ^ 67 is the 
main goal. Therefore, they are already prepared to fight the background source given by the 
radiative decay /i euv"y by using kinematical cuts that suppress undesired signals, although 
reducing the allowed phase space for /j, — > eJ-f as well. Moreover, one has to take into account 
that some experiments, like MEG [S], are only interested in muon decays including photons, 
and therefore they are only sensitive to the process n — > eJ^. For more details, see [35]. 
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Figure 5. Br{fj, — > eJ) as a function of vr for different values of The parameter is taken 
in the range [1,75] TeV, with darker colours indicating a larger value. The MSSM parameters 
have been adjusted in the same way as in figure[2j The points have been selected so that neutrino 
data can be correctly fitted at tree-level. 

5. Summary 

We have studied the phenomenology of a model based on the spontaneous breaking of R-parity. 
In particular, we have concentrated on the decays of the lightest neutralino, showing how we can 
extract interesting information from their correlations with neutrino physics. These correlations 
allow us to test the model at the LHC, confronting the values of special combinations of branching 
ratios with current experimental ranges for the neutrino mixing angles. 

In addition, we have shown the impact of the model on muon decays, which can be 
significantly modified due to the existence of a massless Majoron. Current and future 
experimental facilities might be able to find a deviation from the standard expectations. 

In conclusion, the model is predictive and can be tested at the LHC. 
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